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The magic wavelength for an optical lattice for hydrogen atoms that cancels the lowest order AC 
Stark shift of the 1S-2S transition is calculated to be 513 nm. The magnitude of AC Stark shift 
AE = -1.19 kHz/(10kW/cm=) and the slope dAE/dv = -27.7 Hz/(GHz ■ 10 kW/cm^) at the 
magic wavelength suggests that a stable and narrow linewidth trapping laser is necessary to achieve 
a deep enough optical lattice to confine hydrogen atoms in a way that gives a small enough light 
shift for the precision spectroscopy of the 1S-2S transition. 
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I. INTRODUCTION 

Hydrogen spectroscopy is of fundamental interest to 
physicists, and has contributed to the development of 
quantum mechanics and quantum electrodynamics [l|. 
The spectroscopy has become more and more precise as 
new technologies have developed. In particular, the pre¬ 
cision of the spectroscopy of the 1S-2S transition has im¬ 
proved by ten orders of magnitude in the past century, 
and now the fractional frequency uncertainty is on the 
order of 10“^® ii- This high precision contributes 
to the determination of several fundamental constants, 
such as Rydberg constant and proton radius Q, and is 
also used to set limits on the time variation of the funda¬ 
mental constants @ and the violation of Lorentz boost 
invariance @ . It is also planned to test the CPT theorem 
by comparing the transition frequency of hydrogen and 
anti-hydrogen 

So far, the precision spectroscopy of the hydrogen IS- 
2S transition has been performed with a hydrogen atomic 
beam, because of the difficulty in trapping and cooling 
hydrogen Q . Spectroscopy with an atomic beam cannot 
avoid the uncertainty due to the limited amount of inter¬ 
rogation time and the Doppler effect, and indeed the 2nd 
order Doppler effect is one of the major sources of the 
frequency uncertainty in Ref. Q. The precision spec¬ 
troscopy of other atomic species, on the other hand, is 
typically performed with trapped atoms or ions and takes 
advantage of the long interrogation time and the Lamb- 
Dicke regime confinement, which results in better relative 
uncertainty Eini- For the hydrogen 1S-2S transition, 
one would also expect that spectroscopy with trapped 
atoms would improve the precision. 

To trap neutral atoms for precision spectroscopy, an 
optical lattice formed by a standing wave of laser light is 
typically used. The light-induced AC Stark shift becomes 
a trapping potential for atoms, but since the amount 
of the AC Stark shift is generally different for different 
states, the laser light also induces a frequency shift in op¬ 
tical transitions. At a special wavelength for the trapping 
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light called the magic wavelength, the AC Stark shifts 
for the ground state and an excited state are the same, 
which nearly leads to the cancellation of the energy shift 
of the transition. The idea of the magic wavelength was 
first proposed for the strontium clock transition [l2j | , and 
is now widely used in state of the art optical transition 
atomic clocks la¬ 
in this paper, the magic wavelength for the 1S-2S tran¬ 
sition of hydrogen is calculated, and the possibility of 
trapping hydrogen in an optical lattice of the magic wave¬ 
length is discussed. The AC Stark shift for the hydrogen 
ground state has been widely calculated. However, some 
[l3l - [l^ are for the purpose of high intensity laser appli¬ 
cations, and some others [1^ are calculations in a general 
situation. To the best of my knowledge, there has never 
been a report comparing the AC Stark shift of the ground 
state with that of the 2S state for precision spectroscopy. 

In addition to the 1S-2S spectroscopy of hydrogen, op¬ 
tical trapping is particularly important for anti-hydrogen 
spectroscopy, where an intense atomic beam for spec¬ 
troscopy cannot be generated. The state of the art anti¬ 
hydrogen trap for spectroscopy is a magnetic trap, and 
some measurements were performed for the ground state 
hyperfine transition Q, where the effect of the magnetic 
field is removed by subtracting two frequencies from tran¬ 
sitions between different sub-levels. With a magnetic 
field-free measurement in an optical lattice, the effect 
of the magnetic field is automatically removed, and the 
overall sequence to reduce magnetic field effects becomes 
simpler. 


II. CALCULATION 

To calculate the trapping depth by an optical lattice 
for the ground state and an excited state, I calculate the 
AC polarizability of atoms in those states. Typically, 
such a calculation is performed with relativistic many- 
body perturbation theory [T^, [13) Ell > but for hydrogen, 
simple non-relativistic perturbation theory with the ana¬ 
lytic solution of Schrodinger equation can be used. Some 
of the previous reports on the hydrogen ground state AC 
Stark shift have also used this simple method [l3, H] • 
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The two lowest order energy shifts of a state due to 
the oscillating electric field are given as 

Ai; =-ia(e,w)£:^ - ^7(e,a;)£’^-, (1) 

where a, 7 , e and w are the amplitude of the electric 
field, polarizability, hyperpolarizability, the polarization 
of the light and the light frequency [l^ • The largest con¬ 
tribution to the polarizability is from electric dipole (El) 
transitions, and electric quadrupole (E2) and magnetic 
dipole (Ml) transitions have the second largest contribu¬ 
tion: 

a(e, w) = aEiie,uj) + aMi{e,uj) + 0 ^ 2 ( 6 , w) H- ( 2 ) 

Second-order perturbation theory gives aEiie,uj) for a 
state \n) as 


a_Ei(e, w) 



u)kn\{k\A ■ e|n)|" 


-'kn 


(3) 


The matrix element is separated into the angular com¬ 
ponent and the radial component. The angular com¬ 
ponent involves the Clebsch-Gordan coefficients, and is 
generally magnetic sub-level dependent. In the case of 
the nS-n’P transitions in hydrogen, which are the only 
allowed El transitions from nS states, the transitions be¬ 
tween F — 0 component of a nS state and F = 1 com¬ 
ponent of a n’P state will have polarization independent 
Clebsch-Gordan coefficients. For a transition between a 
^Si /2 F = 0 state and a ^Pi /2 F = I state, the coef¬ 
ficient is l/-\/3, while for a transition between a ^Si /2 
F = 0 state and a ^P 3/2 F = 1 state, the coefficient is 
y^2/3. Since the polarization independent AC Stark shift 
by the trapping light is desired, 1 assume that the spec¬ 
troscopy is performed between F = 0 subl evels of the IS 
and 2S states and regard l/-\/3 and y^2/3 as the angular 
components of the matrix elements respectively. 

The general form of the radial wave function of the 
hydrogen atom is 


where hujkn is the energy difference of the state \k) and 
\n) and d is the operator for dipole moment. For hy¬ 
drogen, the exact energy levels and wave functions can 
be found as the solutions of the Schrodinger equation 
in a nonrelativistic treatment, while relativistic correc¬ 
tion, given by the difference between the solution of the 
Schrodinger equation and the solution of the Dirac equa¬ 
tion, can be regarded as a higher order correction. 1 first 
calculate the lowest order shift, and then estimate the 
corrections due to higher-order terms. 
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where ag is the Bohr radius and 


2Zr 

(5) 
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With Eqs. |4]and[5l the radial components become 
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Combining the angular components, Eqs. [3lEland[71 and setting Z = 1, I obtain 
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where a is the fine structure constant. The summa¬ 
tion over n does not have any simpler analytical form 
and therefore can be calculated numerically with a large 
enough upper limit Umax- 

Figure [ 1 ] shows the calculated AC Stark shift with 
n-max = 100 for visible light. The IS state has an almost 
constant AC Stark shift in this region, while the 2S state 


shift changes considerably. This is because the minimum 
transition energy for the IS state is 10.2 eV, and visible 
fight is far red-detuned for all transitions from the IS 
state. This, in turns, means when the AC Stark shift for 
the IS state changes significantly due to the transition 
between the IS and nP states, the shift for the 2S state 
is more or less constant and positive. Since a negative 
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AC Stark shift is required for the one dimensional optical 
lattice, the visible light region is of interest. 

Figure [T] suggests that there is a point where = 
a ‘^1 around huj ~ 2.4 eV, and Fig. [2] is the magnified 
plot for this region. The crossing point is at 2.4185 eV, 
which is 512.64 nm in the units of wavelength. This is the 
lowest energy magic wavelength for the hydrogen 1S-2S 
transition. Based on the fact that there are shorter wave¬ 
length transitions for the 2S state, there are more magic 
wavelengths for the hydrogen 1S-2S transition, such as 
2.806 eV (441.8 nm) and 2.997 eV (413.7 nm). However, 
2.4185 eV is the best magic wavelength in the sense that 
the slope dAE/dut is the smallest. Thus, I will concen¬ 
trate on the magic wavelength of 512.64 nm. 

III. EFFECTS OF HIGHER ORDER TERMS 



FIG. 1. (Color online) AC Stark shift for the IS (black 
solid line) and the 2S (red dashed line) state by visible light. 
Intensity is lOkW/cm^ and rimax = 100. 


Next, I estimate the contribution by higher order 
terms. The effects of E2 and Ml transitions are calcu¬ 
lated by perturbation theory. The Hamiltonian for the 
E2 transition is 






( 10 ) 


With a linearly polarized plane wave E = 
the polarizability by E2 transition is given as 


a_E2(w) 


2 

h 4c^ 


E 



(11) 

The polarizability of the IS state due to the 3D state is, 
for instance. 
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1S-3D 

E2 


(w) 


3*^ /ao\4 UJ 31 
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( 12 ) 


where the angular component takes a value between 0 
and 4/45 depending on the direction of e and n, and 
4/45 was used to set an upper limit. The polarization of 
the IS state by the 3P state is 


a 


1S-3P 

E\ 


(w) 


3® W 31 

h\z) ’ 


(13) 


and the difference is the factor of , except for the 

numerical prefactor. With visible light and an atom, this 
factor is negligibly small and thus oie 2 is expected to be 
negligible compared to asi- At 512.64 nm, the ratio be¬ 
comes 2.37 X 10“^. Because this factor of w^Oq/c^ is dom¬ 
inant in the ratio the ratio of the po¬ 

larizability due to other excited states c(^e 2 ~^^ / 
is also expected to be small. Thus total polarizability 
due to E2 transitions Q!_e 2 (w) is negligible compared to 

aEi{uj). 

The Hamiltonian for the Ml transitions is 

HMi = f^{L + gS)-B (14) 
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FIG. 2. (Color online) AC Stark shift for the IS (black 
solid line) and the 2S (red dashed line) state around 2.4185 
eV (512.64nm) by lOkW/cm^ light, rimax = 3000. 



and this gives 


aMi(w) 


2 2 ^kn\{k\ {L + gS) ■ e\n)\‘^ 


(15) 


for a linearly polarized plane wave B = 
(i?o/c)ee*^^'®“‘^‘. This is significantly smaller than 
a_Ei(w) due to two factors. The factor is o? smaller 
than c^Oq. The matrix element suggests that there is no 
change in electronic structure, and thus the only allowed 
Ml transitions for hydrogen are hyperfine transitions. 
This results in Wfcn smaller by factor of 10® compared 
to the El transition case. Thus, aMi(w) is negligible 
compared to aEi{uj). 

Hyperpolarizability 7 ( 0 , w) is induced by higher order 
perturbative interactions between an atom and two pho¬ 
tons. Given that the electric field of the light is much 
smaller than the internal field in the atoms, as shown by 
the six or more orders of magnitude smaller AC Stark 
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shift for typical trapping depth than the atomic energy, 
the perturbative expansion is a good approximation. In 
this case, 7(0, cu) is significantly smaller than a(e, w), un¬ 
less there is a two photon resonance. Since neither the 
IS state or the 2S state has a transition of 4.8 eV, no 
two photon transition is expected to give a significant 
contribution to 7(0,0;). 

The relativistic correction is divided into a correction 
of the energy and the correction of the wave function. 
These are expected to be small, as the velocity of the 
electron is ca, and is thus not relativistic. Based on the 
Dirac equation for the hydrogen atom, I can estimate the 
effect with analytic solutions. The energy level correction 
includes a factor of — w^). For example, the rel¬ 

ativistic correction to the IS1/2 state and the 2P3/2 state 
are -43.8 GHz and -2.74 GHz, which gives the change 
of 41.1 GHz in a;2i. This shift in 0721 gives the relative 
change of a;2i/(w|i — of 1.86 x 10“^ at the wavelength 
of 512.64 nm. 

The relativistic correction in the wave function slightly 
changes the matrix element. Defining the difference be¬ 
tween the relativistic wave function and the non- 

relativistic wave function as I obtain 

assuming the correction is small. We calculate 
to estimate As 


Wi^)=YS{e,cj^U2T{l + 2^) 

-iji 


f-V f-V 

Vao/ \ao ) 


^-Zr/ao 


(17) 


where 7 = \/l — the matrix elements are 
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2S 
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IS 
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1.29020^ (18) 

Zj 

128 ciq 

1.29026^ (19) 


The difference is = 5.4 x 

(4'2skl^ir''''') is expected to 
be on the same order of magnitude, and therefore the 
overall relativistic correction to the wave function is 
negligible. 

All other effects on the energy levels, such as the Lamb 
shift, the finite nucleus size effect and the hyperfine split¬ 
ting are around 1 GHz or less, which corresponds to 
a fractional amount of 10“® or less. Thus, the overall 
higher order correction is at most on the order of 10“^, 
and therefore the number 512.64 nm is reliable up to the 
three digit precision. 


TABLE I. Scattering rate by lattice light for relevant transi¬ 
tions 


Transition 

P [S-I] /sat 

[mW/cm^] 

R [s-i] 

IS ^ 2P 

6.26 X 10® 

7256 

7.9 X 10"® 

2S ^ 3P 

2.24 X 10 '^ 

1.65 

1.4 X 10"® 

2S ^ 4P 

9.67 X 10® 

1.75 

1.6 X 10"® 


IV. IMPLEMENTATION OF THE MAGIC 
WAVELENGTH OPTICAL LATTICE 

Figure [2] shows that the AG Stark shift at the magic 
wavelength is -1.19 kHz/(10kW/cm^). This is 50 times 
smaller than the alkali earth-like atoms like strontium or 
ytterbium [l^. In addition, the recoil energy for hydro¬ 
gen is 72 /rK, due to its small mass. These two factors 
require very intense light for hydrogen trapping. In order 
to obtain a trap depth for hydrogen of 300i?r trapping 
depth, which i s ty pical for the state of the art optical lat¬ 
tice clocks [23,[2i|, 3.8x 10® kW/cm^ intensity is required 
for hydrogen trapping. This is barely achievable by fo¬ 
cusing 1 W light injected to an optical cavity of finesse 
3,0007r that has a beam waist diameter of 10 /iin. 

The slope of the polarizability is dAa/di/ = —27.7 
Hz/GHz for lOkW/cm^ light. The dOOAj. lattice gives a 
-10.5 MHz/GHz shift. Gompared to the ytterbium magic 
wavelength trap that gives a slope of 11(1) Hz/GHz for 
SOOF/r lattice [^, this is six orders of magnitude larger. 
In order to suppress the frequency uncertainty of the IS- 
2S transition due to the light shift to IHz or lower, the 
linewidth of the trapping light should be 100 Hz or less, 
and the magic wavelength needs to be determined with 
a similar accuracy. 

Given the high intensity of the lattice, the loss from 
the lattice due to the scattering of the lattice light be¬ 
comes a concern. Table U summarizes the linewidth F, 
the saturation intensity /gat, and R, the scattering rate 
at 512.64 nm for the three closest transitions involving 
the IS and 2S states. The rate is significantly smaller 
than 1 and therefore the loss due to the scattering is 
not a concern. However, it should be noted that the small 
mass of hydrogen and the optical cavity to enhance the 
power might complicate suppression of the heating due 
to the lattice intensity fluctuations [^ . 

Another practical concern is cooling hydrogen to a 
temperature cold enough to trap into the lattice. Hy¬ 
drogen was first trapped in a magnetic trap with buffer 
gas cooling and then evaporatively cooled down to 50 
/iK to achieve a Bose-Einstein condensate [2^. This is 
cold enough for atoms to be loaded into the optical lat¬ 
tice, but this method does not work for anti-hydrogen, 
as the number of atoms that can be trapped is only a 
few in each cycle of the experiment. A cooling scheme 
with the Lyman-a transition was recently proposed , 
and the predicted achievable temperature was 20 mK. It 
would be possible to trap a few anti-hydrogen atoms in 
the optical lattice, but a more sophisticated way to cool 
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anti-hydrogen would be necessary. 

V. CONCLUSION 

The magic wavelength for the hydrogen 1S-2S transi¬ 
tion is calculated using the solutions of the Schrddinger 
equation. The wavelength is estimated to be 513 nm; the 
trap depth is -1.19 kHz per lOkW/cm^ intensity; and the 
slope of the transition frequency shift is -10.5 MHz/GHz 
for a 300 Er lattice depth. These results imply that a 


sophisticated system is required to implement an optical 
lattice for hydrogen, such as power enhancement by an 
optical cavity and a narrow linewidth laser. However, it 
should be possible to trap hydrogen atoms in an optical 
lattice for the 1S-2S transition spectroscopy. 
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